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DISPERSION   OF   EFFLOENT   PLUMES 

FROM  DIFFUSERS    ON  NEAR-SHORE   REGIONS    OF  THE   GREAT  LAKES 

(SURFACE   DILUTIONS) 


SUMMARY 

A   two-dimensional   analytical  model   to   predict   surface   dilutions    from 
a  diffuser   outfall,    in  a    steady   current   regime,    is   discussed. 
Simple  graphical   and   analytical    solutions    are   presented   along  with   a 
fully  worked   out   example.      For  most    cases,    graphical    solutions    are 
recommended    for    simplicity   to   obtain  dilution    ratios    at    any    required 
distance    from   the    source    or   diffuser.      Analytical    solutions   may   be 
used   when    Che   results    are   beyond    the   range    of   the   graphical     , 
solutions. 


DISPERSION   OF   EFFLUENT   PLDMES 

FROM   DIFFDSERS    ON   NEAR-SHORE    REGIONS    OF    THE    GREAT   LAKES 

(SURFACE   DILUTION) 


INTRODUCTION 

After    Che   waste   has    been    initially   diluted    (Handy,    1981)    by    the    jet 
rise,    the  effluent    field   begins    to   drift   with    the   prevailing 
currents    and   spreads    on   the    lake   surface    in   the    shape   of  a   plume  by 
entraining   freshwater   at    the   plume   edges.      This    is    referred    to   as 
surface    dilution. 

Currents 

The  lake  currents  are  generally  variable  with  the  time  (see  example 
for  Lake  Ontario,  Figure  1).   Generally,  the  lake  is  considered 
stagnant  if  th^  current  speed  is  less  than  5  cm/s.   Occassionally, 
the  current  speeds  can  be  as  high  or  higher  than  30  cm/s  due  to  the 
coastal  jetting.   It  is  important  to  use  the  appropriate  value  of 
current  speed  in  any  prediction  model.   If  the  water  movement  data. 
in  Che  region  of  concern,  are  not  available,  current  speeds  and 
directions  should  be  measured  by  drogues  under  different  wind 
cond  irions . 

The  variability  of  the  currents  (see  Figure  1)  in  magnitude  and 
direction  are  responsible  for  the  dispersion  and  dilution  of  a 
pollutant.   In  steady  state  disperion  models,  an  average  speed  is 
generally  used.   Variable  currents  are  utilized  in  more  complex 
time-  variant  models,  such  as  hydrodynamic  and  stochastic  models.   A 
computer  model  (Kuehnel  et  al.  1981),  using  the  variable  currents  or 
the  time  series  of  current  speed  and  direction  over  a  period,  is 
presently  being  tested.   As  this  model  will  use  the  current  meter 
data  (Figure  1)  and  not  just  an  average  current  speed,  it  is 
expected  to  yield  more  reliable  results.   However,  its  use  will  be 
limited  to  the  sites  where  current  meters  have  been  operated  for 
extended  periods. 


Surface   Plume 

Figures    2   and   3    depict    the   behaviour   of   typical   plumes    from  a 
multiport   diffuser.      Due    to   the   jet   rise,    the   plume    from   each  port 
will    travel   a  horizontal   distance   x. ,       from   the   diffuser 
centreline    and   have   individual   plume   widths   w   and   centreline 
dilutions    D.    at    the    surface    (see   Figures    2    and    3).       Computation   of 
X, ,    w  and   D^ ,    in  the   nearfield,    have   been   discussed   by  Hamdy 
(1981). 


The   sewage    field   is    then   further   diluted   due    to   one   or   all    of   the 
following  processes: 

(a)     Midfield   dilution    (D^)    of   individual    surface   plumes,    from 
each   diffuser   port   before    they  merge    into    one    another    (see 
Figure   2). 


(b)      Farfield  dilution    (D-)    of   the    already  merged   plume    (see 
Figures   2    and   3). 


(c)     Apparent   dilution    (D.)    applies    to   the   decay   of 

non- conservative  pollutants,    such   as    coliforms    and   chlorine 
at    all    points    in  the   plume. 

Dilutions    D^  ,    and   D^    and   D,    of   a   conservative    pollutant 
maintains    the    total  mass    of    the    pollutant    but    reduces    the 
concentration   of   the   pollutant.      However,    D      reduces    the   mass    of 
the    non-conservative   pollutant    due    to    the   decay   or  bio-chemical 
reactions;    thereby   reducing    the    concentration   of   such   pollutant. 

DILUTIONS 

Various    centreline    dilutions    are    defined    as    follows;    the    proper   use 
of   svmbols    is    illustrated    in   Figures    2    and    3. 


Nearfield  Dilution,    _      _   C  (1) 


-   2 


where   C^   =   effluent   concentration  of  pollutant  at    the  exit   of  the 
diffuser 
C^    =  plume  centreline   concentration  of  the  pollutant   at   the 
surface,    at   the   interface  NN    :i.e.    interface  of 
nearfield  and   midfield,    Figure  2). 

Midfield  Dilution,  _  C  (2) 

where   C^   =  plume  centreline   concentration  of  the  pollutant  at   the 

interface  FF  of  midfield  and  far- field   (Figure  2  and  3). 

Farfield  Dilution,  _   C  (3) 

3   ■  — 
"2 

where   C^   =  plume  centreline  concentration  of  the   pollutant    in  the 
farfield  at  a  distance  x.    from  the  "interface  FF 
between  mid- field  and  far- field. 
Dj^  =  apparent  dilution  due   to  decay  of  non-conservative 

pollutants.     This  will  be  discussed,    in  more  detail,    in 
a  subsequent   section. 

Total  Dilution,    D^   =   D^    .    D^    .    D         D^  (H) 

at   X  =   X,    +  x^   +  x^    vsee   figure  2)  (5) 

When  there    is   no   midfield  dilution: 

X-    =   0 
D^   =    1.0 
-         C^  =   C^ 

The  midfield  plume  will   develop  only  if  w<f  where      f  =   port 
spacing  and  w,    the   width  of  each   plume.      When  the  plume  width   has 
grown   from  w  to    [     [see  Fig-are   2),    the  midfield  region  ends  and   the 
farfield   region  begins. 
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If  w  <  f ,  there  will  be  no  midfield  region  (see  Figure  3)  as  the 
individual  jet  plumes  would  have  merged  before  reaching  the  lake 
surface. 

Since  the  surface  plume  advection  and  dispersion,  in  the  midfield 
and  farfield  region,  are  governed  by  the  prevailing  currents  and  the 
width  of  the  initial  sewage  field,  one  model  would  be  applicable  to 
both  regions.   The  width  of  the  initial  sewage  field  is  w  for  the 
midfield  and  the  diffuser  length  (b)  for  the  farfield. 

Modelling 

Several  one -dimensional  and  two  dimensional  models  for  predicting 
the  dilution  of  pollutant  concentration  of  the  waste  plume  have  been 
developed.   Foxworthy  et  al  (1966)  developed  a  statistical  model 
requiring  field  measurements.   Brooks  (1960)  approached  the  problem 
analytically  to  solve  the  classical  partial  differential  equations. 
From  the  several  models  available,  the  one  developed  by  Brooks 
(1960)  has  been  chosen  for  this  manuscript  for  its  simplicity  and 
reliability.   The  only  measured  or  estimated  input  required,  for 
this  model,  is  the  average  current  speed.   The  model  is  briefly 
described  below. 

Brooks  Model 

The  following  are  the  main  assumptions  of  this  model. 

(a)  Eddy  diffusivity  (E)  is  a  variable,  and  is  a  function  of 
the  initial  plume  width  (  p )  or  the  diffuser  length  (b). 

(b)  Longitudual  mixing  (along  the  direction  of  current)  and  the 
vertical  mixing  is  negligible. 

(c)  steady  flow  of  discharge  from  the  outfall. 


In  general,    dilution  D,    due   to   the   initial  plume  width    (p)    and   at   a 
distance  x,    from  virtual   source   can  be   estimated   by  Brooks  model. 


D  =  erf^'l.5/r[(l    +  0.57    13  ix/p)  j    '  -   ^1 


(6) 


where  p  =   initial  plume  width 

/3  =  8E/Up  (7) 

E  =   0.01    p  in  cgs  units  (8) 

=  Eddy  Diffusivity 
D  =  mean  speed  of  prevailing  currents 

The   sewage  width   (L)    can  be   computed  as: 

^/p     =     Ti   +  0.67  /S/x\]^-5  (9) 

The  dilutions  as   computed  above   are  for   the   centreline  of  the  plume, 

in  X  directions    (y  =  0).  Thus,    D,  =     D  (already  computed) 

If  dilutions  at  any  point  (x,    y)    in  the  plume   is   required,    it  can  be 
computed  as: 


D,  ,      =      D,  ,.        D_  (10) 

(x,    y)  (x,    o)  5 


where 


D5      =      ex?      [-18  (y/L)   2]  (11) 

Symbols  x,y,    b  and  L  are    illustrated  in  Figures  2   and   3. 

As  analytical  solutions   of  equations    (6),    (9)   and    (10)   may  not   be 
convenient  for   some  users,    the   respective  monographs  are  presented 
-in  Figures   4,    5   and   6.      Dilutions   and   plume  width  at   any  distance 
can  be  obtained   from  graphs   by    first   computing  /3Cx/p)and   for  the 
two-dimensional  dilution,     [  i8(y/L)2J  will  have   to   be 

computed. 
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Equation    (6)    can   always   be   solved   with    the  help   of  Table    1    (see 
Sample   Collections). 

Midfield  Dilution    (D.) 

If  w  <t  ,    the  midfield  region  will    develop.      The    length    of   the 
midfield  region   x_    (see   Figure    2)    can  be   computed   as: 

X2=[l2.5   Uw2/fl-33][cf/.)Û.67    .^J  ^^2) 

where      t    =  port   spacing 

It    is    very    important    to   keep   all    variables    of    (12)    in   cgs   units 
(i.e.    cm  and    sec)    and   obtain  x, 
to   any   other   units,    if   desired. 


(i.e.    cm  and    sec)    and   obtain  x^    in   cm.      x     can    then  be   converted 


For    the  midfield  dilution   D.",    x      should   be  used    in  equations    (6) 
and    (10)    instead   of  x'and  w  instead  of  p   in  equations    (6)    to    (10). 
Note    that   equation    (9)   may   not   be   required    for   the  midfield  region 
because 


Max   (L)    =     {    =  port    spacing 

Intermediate    values    of   plume   width    in   the   midfield   region    are    not    of 
any    significant   use.      If   necessary,    the    width    of    the   plumes    before 
merging  may   be    computed   by    (9)    using    the    appropriate    x   and    p   values. 

Farfield   Dilution    (D.) 

The    intial   plume   width    of   the    farfield    is    generally   equal    to   the 
diffuser    length    (b).      Equations    (6),    (7),    (8),    (9),    (10)    and    (11) 
may   be   applied    to    the    farfield   region   by   replacing 


X   by   X 

p   by   b 
D   by   D, 
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Decay  Dilution    (D, ) 

The    concentration   of  non-conservative    pollutants    are    reduced   bv 
decay   of   their  mass    (by   biochemical    reaction,    mortality   of    indicator 
bacteria,    etc).       Such   decay    from   all    causes    except   dispersion,    may 
be   expressed    by    the    first    order   equation: 

D^   =   exp    (-kt)  (13) 

where   k   =    first-order    decay   rate    to   be    determined   experimentally    or 
may   be   obained    from   literature.      Note    that   equation    (13)    is 
applicable    to   non-conservative   pollutants    (e.g.    Chlorine,    ammonia. 
phenol   and    indicator   bacteria,    etc.)    whose   decay  may   be   approximated 
by    the    first-order   kinetics. 

t   s   X  =   time   of   travel 
u 
Values   of  k   should  be   determined    from   field   studies   where 

possible;    some   tjrpical   values   obtained    from  our   unpublished 
studies   are: 

for   colifonn,    k   =    16x10*   /s 
for   chlorine,    k   =    20xl0~    /s  in   sunmer 

=    6x10      /s  in  winter 

for    a   conservative    pollutant,    k   =   0,    and    therefore    D     =    1.0 

4 

SAMPLE  CALCULATIONS 

Lakeview  Diffuser  Outfall 

A.   Present  case;  Effluent  discharge  rate,  Q   =  1.96  m^ / 

Hamdy  (1981)  have  computed  the  following  for  the  above  example. 


^  =    12   m 

w  =   4.4  m 

x^  =    12  m 

D^  =   0.077 

Since  w<f,    there   is  no   interference  of  individual   jets. Figure  2 
depicts  the  behaviour  of  such   sewage   field.     Thus,    there  will  be 
midfleld  spreading  resulting  in  dilution  D_. 

Consider  the   case  with  a  mean  current  speed,     U  =    10  cm/s. 

(a)   Midfield  Dilution  D^ 


Using  equation    (12),    we   compute   x_ 


Remember   that  all   variables  must   be  in   cgs   (cm  and   s)    units   for  this 
equation;    therefore, 


X2    =  12.5(10)(Ui40)2 


[i^r  - 1 


(1200)^-33 

=   1898(0.95)   =   1807  cm 
=   18.07  m 

Since  we  need  to  determine  dilution  D-  at  x_  =  18  m  for  a  sewage 
width  of  4.4  m,  the  use  of  nomograph  is  not  recommended  because 

/3(x/b)  or  /3(x2/w)  is  less-  than  1.0 

Using  the  more  accurate  analysis  of  equations   (6),    (7),    and   (8)    as 
follows: 


From  (8) 


From   (7) : 


E  =   0.01(440)^'^^    =   33.5    cm^/s 


6     -  8E/Uw  -  8x33-5   =  0.061 
10x440 


X2    =    1800   cm 


From   (6); 


D2   =  erfJll.S/     j"  1    +   0.  67  (0.  061  )/l800\)3   _   i| 


•f     /1. 5/(1. 59   - 


=  erf/(2.5)     =  erf  (1.58) 
=  0.9996  from  Table  1 


(b)    Farfield  Dilution   D- 


b   =   210  m     690    ft    (Hamdy,    1981) 
U   =    10    cm/s      20    fpm    (assumed) 


Assume    chat    the    dilution   D-    is    required    at    X.   =   2100m. 


Frcsn.  (8) 


E     -  ,0,1    (21000)^-^^        -        5792   cm^/s 


From  (7): 

B   '       ^     »  8  X   5792      ^  0.22 

ub  10  X    21000 


For 


X  =    2100  m,/3(|)     =      0.22/2100\ 

V    210/ 


=      2.2 


Using    Figure    4   and    ^  /  x   \   =   2.2,      Dilution   can   be   read 


H-f) 


D^   =   0.38 


The    plume   width    (L)    at      3    =   2100m 
can   be   obtained    by   using   Figure    5. 


^(t) 


^°^     o  /   X   \     =      2.2,  L  =     4.0 

b 

L  =      4b      =     4  X   210 

L  =      840  m 


-  9 


(c)    2-diinensional   dilution 

Dilution   at    a  point  whose    co-ordinates    are    (x,    y) 
i.e.    not   on   centreline 

Let      X     =      2100  m        and      y  =   210  m 


D      (x,    y)      =      D-(x,    o)       .    D- 
For  D^,    compute    18    (Y/L)^      -      Isf-^V     =      1.1 


\840/ 


From   figure   6,   D.      =     0.33 


D^    (2100,    210)  =      0.38   x  0.33 

=      0.125 


Similarly,    dilution    at    any  point   within   the    plume    (i.e.    at    a 
distance  y   from   the   plume   centreline)    can  be   computed. 

(d)   Total  Dilution   at   Plume   Centreline, 

at   x_    »   2100  m  for   conservative   pollutant,   D,    ■   1.0 


°x  =  °1    •    ^2    •    °3    •    °4 


=    (0.077)(0.976)(0.38)(1.0) 
=   0.0286 


where   x  =   x.    +   x.    +   x_ 

=    12   +    18  +   2100 
=   2130  m 


Thus    the    total    dilution   of  a   conservative    pollutant    at    2130  m   from 
the   diffuser    is: 


°2130   =   °-°286 


IC 


(e)   Decay  Dilution   (D^) 


Now  consider  a  non-conservative   pollutant,    e.g.    coliforms 
decay  ."ate   for   coliforms,    k  =    l6x10~  /s 


therefore 


t   _  X     -      ^"1300    =        21300      s 

u      "  10 


kt  =    16x10*^  X  21300 


From   (13) 


=   3.408 


D^   =   exp(-3.408)    =   0.033 


For  coliform,    Dg^^   =   D^    .    D^    .    D..    D^ 


=   (0.077)(0.976)(0.38)(0.033) 
=   0.0009^ 
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3 
B.   Ultimate  case,  Effluent  discharge  rate,  Q-  =  7.0  m  /s 

Hamdy  (1981)  have  computed  the  following  for  the  above  example 

f  =  6.0  m 
X  =  10  m 
D^  =  0.20 

Since  the  nearfield  flow  is  unstable,  individual  jet  plume  would 
merge  before  reaching  the  water  surface.   Therefore,  width  of  the 
plume  at  surface  equals  the  diffuser  length,  b  =  210  m  (see  Figure  3). 


In  this  case,  D.  =  1.0 
X,  =  0 


and    the  midfield  region  does   not   exist, 


The    computations    and   results   for   the   farfield   and  decay   are   similar 

3 
to    the   previous   case   when  Q  »    1.96  m  /s. 


therefore: 


=    10  +   0   +    2100 
=   2110  m 


°2110  =   °1    •    °2    •    °3    •    °4 


»   (0.20)(1.0)(0.38)(1.0) 
=   0.076 


If   the    pollutant   were   a   non-conservative,    e.g.    coliforms.      Decay 
rate    for    coliforms,    k   =    16x10      /s,    time    to    travel    a  distance    of 
x  =   2110  m   is: 


211000   ^  21100      s 
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kt  =   16xl0"^   X   21100   =   3.376 


From    (13): 


For    colifonns: 


D^   =   exp(-3.376)    =  0.034 


°2110   =  °1    •    °2    •    °3    •    °4 


=   (0.20)(1.0)(0.38)(0.034) 
=   0.0026 

CONCLUSIONS    &    RECOMMEND ATI ON S 

The   following    conclusions    can  be    drawn: 


(a)  The  most    important   dilution   is   due   to   initial  mixing   and 
jet    rise    (D   ) 

(b)  The   dilution   due    to   decay    (D,  )    of  non-conservative 

4 

pollutants    is   also   significant. 

(c)  The   surface   dilutions    (D      and  D   )    of   the  mid   and    far 
fields    are    very   small    compared    to    other    dilutions    D      and 
\- 

In   order    to  minimize    the    areas    of   the  mixing    zones,    it    is    important 

to  maximize    the    dilution    (D    )    due    to    initial  mixing    and    jet    rise. 

D^    is    the    only   dilution    that    can  be   controlled   by    changes    in 

design    for    the  multiport   diff users.      Dilutions   D„ ,    D,    and  D, 

2        3  4 

cannot   be    controlled   by   any   design   changes    and   depend   upon  natural 
factors,    such    as,    dispersion   and    advection   characteristics   of    the 
site    and    decay    rates    of    the   pollutant. 


Ir   no    data    for    currents    are    available,    the    dilution   computations    for 
average    current    speeds    of   5,    10,    and    15    cm/s    should   be    undertaken   as 
representative    values.      A   judicious   decision   on    the    appropriate 
results    will    then   have    to  be  made   by    the   user.      Current    data   should 
then  be   obtained    as    soon    as    possible   to    confirm    the   predictions. 
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TABLE  1:   Error  Function  (erf) 


X  erf(x) 

0.0  0.0 

0.1  0.1129 

0.2  0.2227 

0.3  0.3286 

0.4  0.4284 

0.5  0.5205 

0.6  0.6309 

0.7  0.6778 

0.8  0.7421 

0.9  0.7969 

1.0  0.8427 

1.2  0.9103 

1.4  0.9523 
1.6  0.9763 
1.8  0.9891 
2.0  0.9953 

2.5  0.9996 
3.0  0.9998 

1.0000 
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FIGURE    4  :    DILUTION    ALONG    CENTRELINE   OF  A  SEWAGE    FIELD    (Brooks   N.H.,1960) 
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